Solitons and decoherence in left-handed metamaterials 
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We present exact electromagnetic solitary pulses that can be experimentally obtained within 
nonlinear left-handed metamaterials. The effect of pulse decoherence on the modulation instability 
of partially incoherent electromagnetic waves is also investigated. The results may contribute to a 
better understanding of nonlinear electromagnetic pulse propagation in media with negative index 
of refraction. 
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The electromagnetic properties of materials exhibiting 
a negative index of refraction was investigated already 
by Mandelstam in 1945 0, 0. In addition, the idea of 
negative permittivity and permeability and the conse- 
quences of a negative group velocity was pointed out by 
Pafomov and Agranovich and Ginzburg , and then 
discussed by several other authors, e.g. Refs. la la B , l a - 
In recent years such materials have been produced |9ul0|. 
and experimental verification of previous theoretical is- 
sues are currently being studied with great interest (see, 
e -g- | 111 1 1 2l for a review). Quite recently, Lazarides and 
Tsironis [L| presented a model for the propagation of 
electromagnetic waves in nonlinear metamaterials hav- 
ing negative index of refraction (left-handed materials). 
They also presented a study of the possibility of soliton 
formation within such materials. Nonlinear properties 
of left-handed metamaterials have also been analyzed in 
Ref. 0. 

In this Letter, we show that the system presented in 
Ref. 0] also allows for dark/bright soliton pairs, exhibit- 
ing pulse width broadening due to the coupling of the soli- 
tons. Moreover, we analyze the statistical properties of 
the nonlinear system of equations, in particular the effect 
of partial incoherence on the modulational instability of 
partially incoherent electromagnetic waves in left-handed 
metamaterials. The relevance of our investigation to lab- 
oratory (photonic crystals) and astrophysical settings is 
discussed. 

Left-handed materials simultaneously have negative 
permittivity e c s and permeability // e ff, and therefore al- 
low for a negative refractive index. As demonstrated in 
Refs. and 0], it is possible to obtain the expres- 
sions e c ff = e + a|E| 2 and ^ e fi = A* + /5 1 H | 2 , where E and 
H are the electric and magnetic fields respectively, for 
weak fields in left-handed metamaterials. Here e and /i 
are the linear permittivity and permeability, respectively, 
which may take on negative values in metamaterials. We 
note that if we have focusing media (a,/? > 0), the non- 
linear effect may be used to switch between left-handed 
and right-handed character of the material. Moreover, 
for defocusing media (a, < 0), we will have a reinforce- 
ment of the left-handed properties of the material. Thus, 
the nonlinear extension of linear metamaterials may have 
important consequences for future applications. 



Using the expressions for e c s and fi c g, Lazarides and 
Tsironis derived a system of equations for the non- 
linear propagation of an electromagnetic pulse in left- 
handed materials, or 

idflp a + d 2 x il a + {o-aHa? + CT b \^ b \ 2 )^ a = 0, (la) 

id t xp b + d 2 x ij b + (a a \ib a \ 2 + aM 2 )^ b = 0, (lb) 

where the coefficients a a = sgn(a/i) and a b — sgn(/3e) 
can have the values ±1, and the envelope amplitudes 
ipa } b (in Ref. 0] denoted by P and Q) correspond to 
the electric and magnetic field components of the electro- 
magnetic pulse, respectively. Thus, for focusing media, 
if jji, e > 0, we have a a = <r b = 1, while if fi > (<)0, 
e < (>)0, we have a a = +1(-1) and a b = -1(+1). If 
both \x and e are negative, we have a a = a b = — 1. An 
analysis of the system Q was performed in Ref. ^| . It 
was found that Eqs. Q admits stationary solutions of 
the form ipi = ipm exp(— iAujt), with ipio > 0, where the 
phase is determined by Auo = — {(J a \^p a \ 2 + o- b \ip b \ 2 ). We 
now perform a first order perturbative analysis of Eqs. 
(JIJ around the stationary state, writing fa = \$Ii(x,t) + 
fao] exp(— iAu/t), where \^i\ <C fao is a slowly varying 
envelope of the wave modulation. Linearizing Eqs. 
we thus obtain (see, e.g. Ref. |l5|) 

id t ^i+dl%+tpio [a a fa (^ a + #*) + a b fa (V b + = 0, 

(2) 

where the asterisk denotes the complex conjugation. 
Thus, with = it, + ivi, the real and imaginary parts 
of Eq. give 



d t Vi + d 2 x u t + 2fa (a a fa l ou a + o- b fa u b ) - 

d t Ui + d 2 x Vi 



0, 

= 0, 



(3a) 
(3b) 



respectively. With the ansatz Uj,Uj oc exp(iKx — iQt), 
we obtain the dispersion relation 

(fi 2 - K 4 + 2a a iP 2 a0 K 2 ){n 2 -K 4 + 2a b i, 2 w K 2 ) 
-4a a a b ^ 2 b0 K 4 = 0. (4) 

Thus, we have a modulational instability growth rate T = 
— zf2 given by 



T = K 



(5) 
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We note that if a a ,b = — 1, the modulational instabil- 
ity is prohibited, reflecting the well-known fact that one- 
dimensional dark solitons are modulationally stable. 

It was shown in Ref. 0] that Eqs. admit coupled 
bright/bright and dark/dark solitons. In addition, how- 
ever, it is straightforward to show that 

fx — 2Kt\ 

ip a (x, t) = ip a0 tanh exp(iKx - ift a t), (6) 



L 



and 



fx — 2Kt\ 

ipb(x, t) = ipbo sech ( - J cxp(iKx - iil b t), (7) 

are solutions to Eqs. if L~ 2 = (a b ^ — a a ipl ) /2, 
tt a = -K 2 - <J a ip 2 a0 , and fl b = -K 2 - (<r a i/j 2 a0 + <J b xp 2 w )/2. 
These coupled dark/bright solitons are solutions, which 
seem to have been overlooked in Ref. JL3J . The effects of 
the nonlinear coupling between dark and bright solitons, 
such as the width dependence on the signs of <j a ,b and 
the pulse intensities, should be experimentally verifiable. 

The statistical properties of electromagnetic pulse 
propagation in nonlinear metamaterials can be studied 
using the Wigner formalism. Introducing the Wigner dis- 
tribution functions (see, e.g. (l7j l 

Fi(x,t,k) = ±- J dte ik ^*(x + Z/2,t)Mx-ZM, 

where i = a, 6. We then obtain the coupled kinetic (or 
the Wigner- Moyal) equations E HJ 

d t F a +2kd x F a +2{a a \ij a \ 2 +a b \^j b \ 2 )sin (\ d x d k \ F a = 0, 



(9a) 



and 



d t F b + 2kd x F b + 2(a a \iP a \ 2 + a b \ 1 p b \ 2 )sm d x d k j F b = 0, 

(9b) 

from Eqs. (JTJ. Since 

\ipi{x,t)\ 2 = J dkFi(x,t,k), (9c) 

Eqs. (|5J) form a closed set of equations for partially in- 
coherent electromagnetic pulses in left-handed metama- 
terials. 

Next, we consider the first order perturbation of Eqs. 
©. We let Fi(x,t,k) = F i0 (k) + F a (x,t,k) exp(iKx - 
ifii), where Fn <C Fq. Linearizing Eqs. © we obtain the 
dispersion relation 



1 



dk 



a a [F a0 (k + K/2) - F a0 (k - K/2)} 



2K J 1 k - n/2K 

a b [F b0 {k + K/2) - F b0 (k - K/2)} \ 
k-n/2K f 



(10) 



A monochromatic pulse can be analyzed by using the 
ansatz F i0 = IioS(k ~ fc i0 ), where ijo — \ipio\ 2 is the in- 
tensity of the field ipio, in Eq. I|I0|I . On the other hand, if 
there is a stochastically varying phase in the background 
envelope fields ifiio , the corresponding Wigner function is 
given by the Lorentzian distribution 



F l0 (k) 



7T (k- k i0 ) 2 + k? T ' 



(11) 



where kiT is the width of the distribution around fc^o- 
With the distribution (|f I |l . the dispersion relation 1I0|) 
takes the form 



1 = — 



2a a K 2 I a 



[fi - 2(k a0 - ik aT )K] 2 - K* 

2a b K 2 I b 

[fl-2(k b0 -ik bT )K} 2 -K i ' 



(12) 



In many cases, it is reasonable to assume that the statis- 
tical spread is equal and centered around the same value 
for the two envelopes ip a ,W- i- e - fcaO.T = &&o,t = fco,T- 
Then Eq. JTJJ yields 

[Q - 2(k Q - ik T )K] 2 - [K 2 - 2{a a I a + a b I b )} K 2 = 0. 

(13) 

Letting = 2k^K + iT in Eq. I|I3|) . we obtain the growth 
rate T for a t I a + o b I b > K 2 /2. We have 



T = K^2(a a I a + a b I b ) - K 2 - 2k T K, (14) 

which shows the damping effect of the decoherence spread 
kr- Moreover, as noted above, there is no growth in the 
case of <j a ,b = — 1- Equation (|14|> also shows the influence 
of a focusing metamaterial. When /i < I , but e > I, we 
have <y a = — 1 and a b = +1, removing the modulational 
instability altogether when I a > I b . 

To summarize, we have presented nonlinear solu- 
tions for electromagnetic fields in materials with nega- 
tive index of refraction. In addition, we have also dis- 
cussed the modulational instability of partially incoher- 
ent electromagnetic waves in such materials by employ- 
ing the Wigner-Moyal equation. Choosing the equilib- 
rium Wigner function in the form of a Lorentzian distri- 
bution, corresponding to a random phase, we have de- 
duced a nonlinear dispersion, which exhibits an oscilla- 
tory modulational instability whose growth rate is given 
by Eq. (|I4|I . The present results should be useful for 
understanding the nonlinear propagation of finite ampli- 
tude electromagnetic pulses in left-handed metamaterials 
in laboratory and astrophysical (e.g. rotating black holes 
plf) environments. For example, multiple black holes 
and other massive objects in outer space can make the 
light beams bend in unexpected and unpredictable ways. 
The unexpected light bending can be attributed to the 
nonlinear propagation of intense electromagnetic pulses 
in a medium with negative index of refraction. 
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